Disentangling density and temperature effects 
in the viscous slowing down of glassforming liquids 
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We present a consistent picture of the respective role of density (p) and temperature (T) in 
the viscous slowing down of glassforming liquids and polymers. Specifically, based in part upon 
a new analysis of simulation and experimental data on liquid ortho-terphenyl, we conclude that a 
zeroth-order description of the approach to the glass transition should be formulated in terms of a 
temperature-driven super- Arrhenius activated behavior rather than a density-driven congestion or 
jamming phenomenon. The density plays a role at a quantitative level, but its effect on the viscosity 
and the a-relaxation time can be simply described via a single parameter, an effective interaction 
energy that is characteristic of the high-T liquid regime; as a result, p does not affect the "fragility" 
of the glassforming system. 



I. INTRODUCTION 

Why should one care about the respective role of 
the various external control parameters, temperature T, 
pressure P, density p (or volume V) in the viscous slow- 
ing down of glass forming liquids and polymers? If the 
glass transition at T g were a bona fide critical point with 
a diverging correlation length, the question of how it is 
approached in the P — T or p — T phase diagram would 
be interesting but not crucial to determine the main fea- 
tures of the critical slowing down. However, there is 
no directly observed critical point associated with the 
glass transition of liquids and polymers. For the canoni- 
cal fragile glassformer ortho-terphenyl, the often invoked 
"ideal glass transition temperature" obtained by extrap- 
olating the Vogel-Fulcher-Tammann expression for the 
viscosity and the a (or primary) relaxation time to infin- 
ity or extrapolating the excess entropy to zero, is some 
40 K below T 9 [1 f 



critical points" at temperatures above j. g , 



Various theories predict "avoided 
T g , such as the 

mode-coupling singularity at T c [3j or the frustration- 
limited-domain avoided transition at T* Q| , or "unreach- 
able critical points" at temperatures below T g , such as 
the ideal glass transition at T a [H El 01 j but whatever the 
approach, the putative correlation lengths are expected 
to remain of modest size, no more than 5 to 10 molecular 
diameters 00. The question of which control parameter 
dominates the slowing down of the relaxation may then 
be relevant to better understand the underlying physics. 

Actually, most zeroth-order theoretical descriptions of 
the slow dynamics of glass forming liquids are formulated 
either at constant temperature or at constant density. 
Two such limiting approaches are: 

1 ) the congestion, or jamming, picture of the glass tran- 
sition, in which slowing down results from the drainage 
of free volume as density increases ■ The paradigm 
for this approach is the hard-sphere fluid; 



2) a description in terms of thermally activat ed p rocesses 
on a constant-density energy landscape [T(J, EH ■ hi this 
approach the landscape does not change as one cools 
the system , and the super- Arrhenius T-dependence of 
the a-relaxation time and the viscosity is attributed to 
changes in the minima and barriers encountered in the 
exploration of the landscape as one lowers the tempera- 
ture. 

Of course, such zeroth-order descriptions can always 
be improved to account for some additional dependence 
on the "secondary" control variable: temperature can 
be included in the hard-sphere picture by allowing the 
diameter of the spheres to change with temperature, and 
the energy landscape can be studied for various densities. 
But the mechanism responsible for the viscous slowing 
down is determined predominantly by the density in the 
former case and the temperature in the latter one. 

In this article, we extend our previous work [Hill to 
develop a consistent picture of the roles of density and 
temperature. To do so, we consider recent studies made 
by others as well as new simulation and experimental 
data. The main conclusion is that the appropriate zeroth- 
order description of the slowing down of relaxation as 
one approaches the glass transition is not that of density 
(or volume)-driven congestion or jamming, but that of a 
temperature-driven super- Arrhenius activated behavior. 
Density plays a role at a quantitative level: there is a sig- 
nificant density dependence of the a-relaxation time and 
the viscosity at fixed temperature. However, the density 
dependence can be described in terms of a single interac- 
tion energy that is characteristic of the high-T Arrhenius 



1 For a given number of particles N and a given volume V, hence 
a given density p = N/V, the potential energy landscape, i.e., 
the hypersurface of potential energy plotted as a function of the 
3N configurational coordinates, is independent of temperature. 



liquid. As a result, the temperature dependence of the 
a-relaxation time and the viscosity at different densities 
collapse onto a single, material-dependent curve; most 
importantly, this implies that density does not affect the 
"fragility" of a glass former, i.e., the degree to which it 
exhibits super- Arrhcnius temperature dependence. 



II. ASSESSING THE INFLUENCE OF DENSITY 
AND TEMPERATURE 

Whereas physical pictures and theories are conve- 
niently formulated in terms of temperature and density, 
experiments are usually carried out by using tempera- 
ture and pressure as control variables. Standard studies 
of the slowing down of relaxation leading to the glass 
transition are performed along isobars, typically at at- 
mospheric pressure, by cooling the liquid or the polymer. 
The resulting T-dependence of, say, the a-relaxation time 
r a includes both an intrinsic temperature effect and a 
density effect resulting from the fact that cooling a sys- 
tem at constant pressure usually increases its density. 
The temperature typically decreases by a factor of 2 be- 
tween the normal liquid and the glass transition, whereas 
in the same range the density increases typically by only 
10%; but of course those numbers are not meaningful by 
themselves, the important physical information being the 
sensitivity of the slowing down to changes of T and p. 

In order to disentangle the effects of temperature and 
density at constant pressure and provide a quantitative, 
model-free assessment of their relative importance, we 
suggested some years ago to use as a criterion the 
ratio of two coefficients of expansivity, \a T \/ap, where 
ap = —p~ l (dp/dT)p is the usual isobaric expansivity 
and a T = —p^ 1 (dp/dT) T is the (unusual) isochronic ex- 
pansivity for a constant a-relaxation time or a constant 
viscosity (77). Indeed, along an isobar the variation with 
T of ln(r a ) where r a is expressed, say, in seconds (or 
ln(ry) with r\ expressed, say, in Poise) can be decomposed 
as 



91n(r a ) 



dT 



91n(r a ) 



dT 



<91n(r Q ) 



dp 



dp 
,8T 



(1) 



where the first term of the right-hand side is the intrinsic 
effect of temperature at constant density and the second 
one is the effect of density at constant temperature. By 
using the above mentioned coefficients of expansivity, it 
is easy to show that the former term is expressed in terms 
of a T and the latter in terms of ap, so that 



91n(r a ) 



dT 



91n(r Q ) 



d\n(p) 



(-a T + a P ), 



(2) 



and the ratio \a T \/ap measures the relative importance 
of temperature versus density contributions at a given 
thermodynamic point (recall that in general ap > and 
a T < 0): if \a T \/ap > 1, temperature is a more impor- 
tant controlling factor of the viscous slowing down than 
density. 
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FIG. 1: Ratio |a r |/ap for various liquids and polymers eval- 
uated at the glass transition and at higher temperature (at 
atmospheric pressure): filled symbols for T g (see text), dotted 
open symbols for r Q ~ 10 -5 s or r a ~ 10 -1 s , open symbols 
for r a ~ 10 -9 s or -q ~ 10 4 cPoise. The two ellipses indicate 
points corresponding to liquids of the same family with a very 
similar chemical structure. The point for DBP at Tg should 
be taken with caution because of the inaccuracy of the equa- 
tion of state. The rank of the liquids (see the table) within 
the three classes, van der Waals molecular liquids, polymers, 
and H-bonded molecular liquids, is arbitrarily chosen. Ref- 
erences for the table: (a) Ref. [12| and references therein, 
(b) Refs. [23, [2land O, (c)Ref. [HJ, (d) Ref. Je) 
Ref. 0, (f) ReFIH, (g) Ref. H, (h) Ref. (i) Ref. 



Note that an equivalent characterization of the respec- 
tive weights of temperature and density, although not 
expressed in terms of thermodynamic coefficients of ex- 
pansivity, was proposed by Naoki et al. |14| . These au- 
thors suggested use of the ratio Ey/Hp, where Ey = 
(dln(T a )/d(l/T)) p and H P = (dln(T a )/d(l/T)) P . It is 
easy to show that 



Ey 

Jlp~ 



ap 



(3) 



so that, when \a T \/ap > 1, Ey /Hp > 1/2. 

From available experimental data on viscosity and a- 
relaxation time we have calculated (or collected) the ratio 
|a T |/ap for a variety of glassforming liquids and poly- 
mers interacting through van der Waals forces or hydro- 
gen bonds. Strong, network-forming systems are not con- 



sidered here 2 . We show in Fig. ^ the ratios calculated 
at atmospheric pressure (i) at T g (either the calorimetric 
T g corresponding to an isochrone with r a ~ 1000 s or the 
dielectric T g with r a ~ 1 — 100 s) and (ii), when possi- 
ble, for a smaller a- relaxation time or viscosity (i.e., a 
higher temperature). It is quite noticeable that the ratio 
at T g is significantly larger than 1, with the exception 
of a family of organic liquids with very similar chemical 
structure (BMPC, BMMPC) for which it is slightly less 
than 1 0. The largest value is for the f/-bonded liquid 
glycerol; van der Waals molecular liquids tend to have 
smaller ratios than polymers and H-bonded liquids, but 
the correlation between the value of the ratio and the 
nature of the glassforming system is not very strong. 

A few additional comments are worth making. First, 
there may be significant errors in determining a T and 
ap: high- precision thermodynamic and dynamic data are 
needed to obtain such second-order quantities. Secondly, 
and maybe somehow contrary to the expectation that 
density becomes more important as pressure increases, 
the ratio \a T \/ap is found to increase with pressure (ex- 
cept for glycerol), at least up to a few kbars. Finally, the 
ratio decreases as the temperature increases or, equiva- 
lently, as r Q (or -q) decreases. It is indeed expected that 
in the "normal" , high-T liquid density and temperature 
play more comparable roles, at least for weakly interact- 
ing systems. However, this is rather at odds with the con- 
gestion/jamming picture of the glass transition in which 
one would expect density to become predominant as one 
gets close to T g . 

To summarize the results of this section: overall, tem- 
perature dominates over density in driving the viscous 
slowing down of glassforming liquids and polymers; for 
a quantitative description one cannot bluntly neglect 
the effect of density; the role of temperature becomes 
more important in the super- Arrhcnius regime as one ap- 
proaches the glass transition. 



III. SCALING OUT THE DENSITY 
DEPENDENCE 

The conclusions of the preceding section, together with 
previous work 0, ^| , lead us to describe the a relax- 
ation in supercooled liquids and in polymers as thermally 
activated, i.e., for practical purposes describable by the 
following T-dependence: 

T a (p,T)=T™(p)exp[E(p,T)/T], (4) 

where the effective activation energy E(p, T) is at least 
equal to several times the thermal energy T (we work in 



2 It is likely that for such strongly interacting glassformers temper- 
ature is the dominant driving parameter of the viscous slowing 
down so long as pressure does not change the local structure of 
the liquid. 
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FIG. 2: Simulation data on a model ortho-terphenyl: 
\n{T/D(p, T)), where D(p,T)) is the translational diffusion 
constant, versus 1/T for several values of p (from top to 
bottom, p = 1.135, 1.108, 1.083, 1.060, 1.036 g/cm 3 ). The 
straight line is an Arrhenius fit to the high-temperature por- 
tion of the data (T < T*(p)), shown here for p — 1.036 g/cm 3 . 

units such that ks = 1). A zeroth-order physical picture 
of the viscous slowing down of fragile glass formers and of 
their crossover from Arrhenius to super- Arrhenius behav- 
ior should thus focus on the role of temperature, the effect 
of density being incorporated in quantities characterizing 
the high-T, "normal" liquid phase. As we suggested in 
Ref. [T3I, the simplest description is that density only 
enters via a characteristic interaction energy that one 
can associate with the energy barrier to a-relaxation and 
viscous flow in the high-T Arrhenius- like regime, Eoo^p); 
this energy sets the scale for both the effective activa- 
tion energy E(p, T) and the temperature. Accordingly, 
all isochoric data should collapse onto a master curve 

= * (—) (5) 

E^p) \Eooip))' {> 

where $(C) is a monotonously decreasing function of £ 
that is equal to 1 for £ larger than some crossover value 
£*. If, moreover, the p-dependence of the prefactor 
in Eq. (0J is negligible, one can also obtain a mastercurve 
for the a-relaxation time itself, 

h(Ta(p ,D) = In(0 + ^*(^), W 

where T a (p, T)and are both expressed, say, in sec- 
onds. It is clear from Eqs. JSJ and © that the effect of 
density could be quantitatively significant without qual- 
itatively determining the slowing down of the relaxation 
as one ap pro aches the glass transition. 

In Ref. [13| we gave two examples of data collapse onto 
the mastercurve in Eq. |5jl (for the simulation data on the 
binary Lennard- Jones mixture [lfij and for experimental 
data on glycerol 0, 0|). It is worth stressing that, 
despite decades of pressure studies of glassforming liquids 
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FIG. 3: High-T Arrhenius activation energy Eoo versus den- 
sity for the data shown in Fig. [5] 



and polymers, the available data base of isochoric data 
is scant. In what follows we present new simulation data 
on a model ortho-tcrphenyl and, building upon the work 
of Refs. ^3 and [2(| we also analyze experimental data 
of ortho-tcrphenyl. 



IV. ANALYSIS OF SIMULATION AND 
EXPERIMENTAL DATA ON 
ORTHO-TERPHENYL 

We have considered the simulation data of Ref. |2l| 
for the now widely studied Lewis- Wahnstrom model of 
ortho-terphenyl 22], a model of rigid three-site molecules 
with intermolecular site-site Lennard-Jones interactions. 
In view of the analysis below, the data base of state 
points has been further extended in the high-temperature 
region. Simulation details are given in Ref. [21|. The 
transnational diffusion constant D(p, T) has been consid- 
ered for 5 different densities in the range 1.036 to 1.135 
g/cm 3 (representing a 10% overall change in density) 3 . 
To obtain data more akin to viscosity or relaxation time, 
we have considered T/D(p,T): see Fig. [5] From the 
Arrhenius plot of Fig. we have determined for each 
curve a high-T activation energy E oa (p) that describes 
the data above some crossover temperature T*(p), below 
which the curves markedly deviate from an Arrhenius de- 
pendence. The values of E^p) are shown in Fig. EI As 
can be seen from Fig.0] the data taken along the 5 differ- 
ent isochores collapse onto a single curve when energies 
and temperatures are rescaled by the parameter Eoa^p); 
the prefactor being independent of density, a master 
curve is obtained by using either Eq. (J5J (Fig. 0] Top) or 




3 Note that, in the range of temperature accessible to our com- 
puter simulations, there is no significant decoupling between the 
translational diffusion constant and the other relaxation times. 



FIG. 4: Mastercurves for the simulation data on model ortho- 
terphenyl at 5 different densities: Top: E(p, T) / E x (p) versus 
T/E x {p); Bottom: ln(T/D(p,T)) versus E^^/T. Symbols 
as in Fig. [5] 



Eq. © (Fig. H Bottom). 

Computer simulation offers a convenient means to 
study isochoric data directly, but the range of temper- 
ature and of relaxation time attainable is quite limited. 
Ideally, isochoric experimental data would provide a more 
stringent test for Eqs. JSJ or @). However, even for the 
well-studied liquid ortho-terphenyl, obtaining isochoric 
data for the viscosity or the a-relaxation time requires 
significant manipulation of the original data (collected 
along isobars and isotherms). We follow here the work 
of others [2(j who gave evidence that the a-relaxation 
time of ortho-terphenyl for a large range of pressure and 
temperature, hence of density and temperature, can be 
placed to a good accuracy on a single master curve when 
it is plotted as a function of pT -1 / 4 : see, e.g., Fig. 11 
of Ref. |2(j ■ We have repeated a similar analysis for the 
viscosity, extending the range to higher temperatures by 
using the data of Schug et al. j^. The resulting mas- 
tercurve, ln(?7) versus p 4 /T is shown in Fig. [S] (the inset 
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FIG. 5: Mastercurve experimental data on the viscosity of 
liquid ortho-terphenyl: log 10 (»7) versus p 4 /T. The data at 
atmospheric pressure are taken from Refs. |33l |34| and those 
at higher pressures from Ref. |23j| (only pressures up to 5 
kbars are considered). The straight line is the best Arrhenius 
fit to the high-T/low-p portion of the data. 



shows a similar plot for the a-relaxation time) . Data col- 
lapse is very good. One should however keep in mind a 
few words of caution. First, there is a significant discrep- 
ancy between cv-relaxation times obtained by different 
techniques. Secondly, we have included data only up to 
5 kbars (similarly in Ref. [2(| only data up to 2 kbars 
were considered); for higher pressures, the viscosity data 
start to deviate from the mastercurve, which may be due 
to a much more uncertain equation of state, to experi- 
mental errors, and/or to the fact that the pT -1 / 4 scaling 
somewhat breaks down at very high pressures. 

What is the meaning of the pT -1 / 4 scaling'? In 
Ref. ^ was suggested that this scaling is reminiscent 
of that occurring in the fluid of soft spheres interacting 
via an r -12 potential. Indeed for particles interacting 
via a power-law r~ n pair potential, the p and T depen- 
dences of all excess thermodynamic properties 4 enter 
only through a single parameter V = pT~ 3 / n [24|; the re- 
laxation time is then essentially a function of T whereas 
the equation of state takes the form PT~( 1+3 /™) = 
F(pT~ 3 / n ) . It would be, however, somewhat surpris- 
ing if a molecular liquid such as ortho-terphenyl behaved 
like a fluid of soft repulsive spheres (with n = 12 as in a 
Lennard- Jones potential) . Actually, from the data in the 
literature it is easy to check that the equation of state 
does not obey the predicted form PT -5 / 4 = T{pT~ x ^). 

We suggest instead that the empirical scaling found for 
the viscosity and the cv-relaxation time of ortho-terphenyl 



4 An excess thermodynamic quantity is defined here as the contri- 
bution due to the interactions that comes on top of an ideal-gas 
result. 
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FIG. 6: High-T Arrhenius activation energy Eoo versus den- 
sity from the viscosity data of ortho-terphenyl. The dot- 
dashed curve is the p 4 -dependence derived from Fig. |^| the 
dotted line is a linear fit to the p 4 -curve over the relevant 
range of liquid densities (shaded area) . 



is one more manifestation of the mastercurve discussed in 
the preceding section: the data collapse shown in Fig. [S] 
is a special case of Eq. © with E 00 (p) oc p 4 . One can 
indeed see from Fig. that the small- (1/C) portion of 
the curve, where 1/C = p 4 /T, is linear to a good approx- 
imation, which corresponds to an Arrhenius behavior, 
]n(r)(p,T)) = 111(7700) + E OQ (p)/T 1 with E oc (p) propor- 
tional to p 4 . In Fig. HO we have plotted E OD (p) for the 
range of densities actually studied experimentally: it is 
obvious from the figure that for this range of densities, a 
linear fit to (p) would do just as well as the power law 
p 4 . This latter may then just be a coincidence without 
any substantial physical meaning. On the other hand, 
the observed data collapse supports our claim that the 
influence of the density on the slow dynamics of glass- 
formers can be reduced to a single effective interaction 
energy E^ip). 



V. COMPARISON BETWEEN SYSTEMS: 
RELATIVE "FRAGILITIES" 

The scaling function $(£) appearing in the master- 
curves given in Eqs. © and © is species-specific, as 
are the p-dependence and the magnitude of the effective 
interaction energy E oa (p). This is illustrated in Figs. 
and |H1 for the two systems studied in the present work, 
computer simulated and real ortho-terphenyl, as well as 
for the two systems considered in a previous study, the 
binary Lennard-Jones mixture and (real) glycerol. In 
Fig. we have plotted Eoo(p) / Eoo(p re f) versus p/pref, 
where p re j is a liquid density taken as a reference (e.g., 
the density at the boiling point when it is known), for 
the four glassformers. One observes that besides the fact 
that over the range considered (a variation of about 10%) 
the p-dependence is essentially linear, the four curves are 
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liquid density, for the glassformers: binary Lennard- Jones 
mixture, glycerol, Lewis-Wahnstrom model of ortho-terphenyl 
(see Fig. [3J , and ortho-terphenyl (see Fig. |SJ . 



all different 5 . 

In Fig. [S] we show the curves E(p,T)/ E OD (p) versus 
ln(T/E 00 (p)) and T/T*(p), where T*(p) = C*£co(p) is 
the crossover temperature below which a marked super- 
Arrhenius T-dependence is observed; in terms of master- 
curves (see Eq. ©), this amounts to plotting $(£) versus 
C/C* where $(C) ~ 1 for C < C*; C* being a species- 
dependent constant. The differences observed for C > C* 
among the four systems reflect the different degrees of 
super- Arrhenius behavior, i.e., the different fragilities: 
one finds, as expected, that (real) ortho-terphenyl and 
glycerol are significantly more fragile than the two com- 
puter simulated systems (the highest points in activation 
energy for glycerol should be taken with caution 13] ). 

It is worth pointing out that a consequence of the exis- 
tence of the mastercurves is that fragility is independent 
of density. This applies no matter what particular defi- 
nition of (dynamic) fragility [2f| one chooses. Consider, 
for instance, the commonly used steepness index, 



dHr a (p,T)) 



d(T g (p)/T) 



(7) 



Tg{p) 



evaluated along an isochoric path. Applying Eq. © 
and defining as usual T g (p) as the temperature at which 
r a (p, T), expressed in seconds, is equal to a given value, 
h\(r a (p, T)) = IC, one finds that 



E^p) d\n(r a (p,T)) 



T g (p) d(T g {p)/T) 



Tg(j>) 



(8) 



This, as well as the results in Fig. [S] casts some doubts on the 
ability of the Lewis-Wahnstrom model to accurately represent 
liquid ortho-terphenyl. 
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FIG. 8: E(p,T)/Eoo(p) versus (Top) ln{T/E x (p)) and (Bot- 
tom) T/T*(p) for the 4 glassformers indicated in the text 
and in Fig. Q T*(p) = (^E 00 (p) is the crossover tempera- 
ture above which the T-dependence of n or r a is essentially 
Arrhenius-like. is species-specific: = 0.42 for the binary 
Lennard- Jones mixture; = 0.075 for glycerol; = 0.53 for 
the model ortho-terphenyl; = 17 for real ortho-terphenyl. 



where ( g = T g {p)/E 00 (p) and d($(C)/C)/dC is nega- 
tive. According to the above definition of T g (p) and 
to Eq. ©, Cg is simply the solution of the equation 
(^KCsVCg) = ~ l n (C)i an d is thus independent of 
p. As a result, whatever the value of JC chosen, the 
steepness index m is independent of the density. This 
remains true if the prefactor has a small residual p- 
dependence. The opposite conclusion reached in Ref. |2|j 
from data on the binary Lennard- Jones mixture proba- 
bly arises from the large intrinsic inaccuracy associated 
with fitting a 3-parameter Vogel-Fulcher-Tammann for- 
mula to a non-fragile system (see Fig. [8J over a limited 
time-scale domain when the putative T Q is at much lower 



7 



temperature than accessible by computer simulation. 

The fact that fragility is independent of density does 
not imply, however, that it is independent of pressure. 
Fragility measured along isobaric paths involves addi- 
tional effects, associated with the variation with p of 
Eoo and of the isobaric coefficient of expansivity ap. 
As a result, fragility may increase, decrease, or be non- 
monotonic as a function of pressure. We will discuss this 
point in more detail in a forthcoming publication. 

VI. CONCLUSION 

Extending our earlier work we have clarified the re- 
spective role of density and temperature in the viscous 
slowing down of glassforming liquids and polymers. On 
the basis of available experimental data, we have argued 
that this slowing down is best described at zeroth- order in 
terms of thermally activated processes whose character- 
istic super- Arrhenius behavior is not primarily driven by 
a congestion or jamming phenomenon resulting, for in- 
stance, from a drainage of free volume as one approaches 
the glass transition, but should rather be explained as an 
intrinsic temperature effect operating at constant den- 
sity. Although density effects cannot be neglected at a 
quantitative level, we have brought additional evidence 
showing that these effects can be described via a single 



effective interaction energy Eoo (p) that characterizes the 
relaxation of the liquid in its high-temperature range. In 
the range of accessible liquid densities, the variation of 
Eoo with p is essentially linear and featureless. 

As a result of the collapse of the viscosity and a- 
relaxation time data obtained at various densities onto a 
mastercurve, we have shown that the (dynamic) fragility 
of glassforming liquids, i.e., their degree of departure 
from the high-T Arrhenius behavior, is independent of 
density. This, in turn, supports the above conclusion that 
a physical explanation of the super- Arrhenius behavior of 
the viscosity and a-relaxation time as one approaches the 
glass transition should be found in the intrinsic effect of 
temperature. Although we feel that we have provided 
strong evidence to support our conclusions, is should 
again be stressed that there is a need for more extensive 
sets of experimental data and accurate equations of state 
to test on a wider scale the proposed scaling form of the 
density dependence of the viscosity and the a-relaxation 
time. 
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